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We study the equilibrium properties of a monomer-monomer A + B — > reaction on a two- 
dimensional substrate containing randomly placed catalytic bonds. Interacting A and B species 
undergo continuous exchanges with particle reservoirs and react instantaneously as soon as a pair of 
unlike particles is connected by a catalytic bond. For the case of annealed disorder in the placement 
of catalytic bonds the model is mapped onto a general spin S = 1 model and solved exactly for the 
pressure in a particular case. At equal activities of the two species a second order phase transition 
is revealed. 
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Catalytically activated reactions (CARs) involve par- 
ticles which react only in the presence of another agent 
- a catalyst, and remain chemically inactive otherwise. 
These processes are widespread in nature and used in a 
variety of technological and industrial applications . 

Following the seminal work of Ziff, Gulari, and Bar- 
shad (ZGB) |2| on the so-called "monomer-dimer" model, 
as well as subsequent studies of a simpler "monomer- 
monomer" reaction model [30, there has been consid- 
erable progress in the understanding of CARs proper- 
ties. The ZGB model, which describes, in particular, 
the important process of CO oxidation on a catalytic 
surface, revealed several remarkable features |2|. On a 
two-dimensional (2D) substrate, upon lowering the CO 
adsorption rate the system undergoes a first-order phase 
transition from a CO saturated inactive phase into a re- 
active steady-state, followed by a continuous transition 
into an 02-saturated inactive phase, which belongs to the 
same universality class as directed percolation and the 
Reggeon field theory [j| . The monomer-monomer model 
exhibits a first-order transition from a phase saturated 
with one species to one saturated with the other; allowing 
desorption of one species leads to a continuous transition 
which also belongs to the directed percolation universal- 
ity class 4] . For these two models, different aspects of the 
dynamics of the adsorbed phase have been investigated 
SSIlilES, confirming an essentially collective be- 
havior. In contrast, the equilibrium properties of CARs 
are much less studied and the understanding of the equi- 
librium state remains rather limited. Recently, an exact 
solution for a one-dimensional monomer- monomer model 
has been presented 0, but for the physically important 
2D substrates no exact solutions are known as yet. 

Although realistic substrates are typically disordered 
and the actual catalyst is an assembly of mobile or lo- 
calized catalytic sites or islands only a few studies 
have addressed the behavior of CARs on disordered sub- 



strates focusing on adsorption/desorption processes but 
not on the coverage of a catalyst decorated substrate [j| . 
With the exception of few exactly solvable ID models 
of A + A — > reactions ^3 an d a Smoluchowski-type 
analysis of ci-dimensional CARs , for random spatial 
distributions of the catalyst only empirical approaches 
based on phenomenological generalizations of the mean- 
field "law of mass action" have been proposed so far pj. 
Consequently, an exact analytical solution of (albeit ide- 
alized) models involving a 2D random catalytic substrate 
is very desirable since it provides valuable insight into the 
effects of disorder on the CARs properties. 

In the following we present such an exactly solvable 
model of a monomer-monomer A + B — ► reaction on 
a 2D inhomogeneous, catalytic substrate and study the 
equilibrium properties of the two-species adsorbate. The 
substrate contains randomly placed catalytic bonds of 
mean density q which connect neighboring adsorption 
sites. The interacting A and B (monomer) species un- 
dergo continuous exchanges with corresponding adjacent 
gaseous reservoirs. A reaction A + B — > takes place 
instantaneously if A and B particles occupy adsorption 
sites connected by a catalytic bond. We find that for 
the case of annealed disorder in the placement of the 
catalytic bonds the reaction model under study can be 
mapped onto the general spin 5 = 1 (GS1) model [l2j |. 
This allows us to exploit the large number of results ob- 
tained for the GS1 model 01 m order to elucidate the 
equilibrium properties of the monomer-monomer reac- 
tion on random catalytic substrates . Here we concen- 
trate on a particular case in which the model reduces to 
an exactly solvable Blume-Emery-Griffiths (BEG) model 
Il5| and derive an exact expression for the disorder- 
averaged equilibrium pressure of the two-species adsor- 
bate. We show that at equal partial vapor pressures of 
the A and B species this system exhibits a second-order 
phase transition which reflects a spontaneous symmetry 



Typeset by REVTeX 



2 



breaking with large fluctuations and progressive coverage 
of the entire substrate by either one of the species. 

We consider a 2D regular lattice of N adsorption sites 
(Fig. [T|>. which is in contact with the mixed vapor phase 
of A and B particles. The A and B particles can adsorb 
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FIG. 1: 2D lattice of adsorption sites (small grey circles) 
in contact with a mixed vapor phase. Black and white cir- 
cles denote A and B particles, respectively. The solid lines 
denote "bonds" with catalytic properties. (R): configuration 
in which an instantaneous reaction takes place upon 
which the reactants leave the system. (NR): NN pair of A 
and B which do not react since the sites are not connected 
by a catalytic bond. 

onto vacant sites, and can desorb back to the reservoir. 
They are characterized by chemical potentials [ia and /ib 
which are maintained at constant values and measured 
relative to the binding energy of an occupied site, so that 
Mas > corresponds to a preference for adsorption. The 
A and B particles have hard cores prohibiting double 
occupancy of sites and nearest-neighbor (NN) attractive 
A — A, B — B, and A — B interactions of strengths J a, 
Jb, and Jab, respectively. The occupation of the i-th 
site is described by a pair Cj of Boolean variables and 
rrij such that 



c l = [ni,mi) 



(1.0) , site i occupied by A, 
(0, 1), site i occupied by B, 
(0, 0), site i empty, 

(1.1) , excluded by hard cores. 



(1) 



We assign to some of the lattice bonds (solid lines 
in Fig. ^) "catalytic" properties such that if an A and 
a B particle occupy simultaneously NN sites connected 
by such a catalytic bond, they instantaneously react by 
forming a product (0) which immediately desorbs and 
leaves the system; A and B particles occupying NN sites 
not connected by a catalytic bond harmlessly coexist. 
The "catalytic" character of the lattice bonds is described 
by Boolean variables C<ij> > where < ij > denotes a pair 
of neighboring sites i and j, 



c 



1, < i j > is a catalytic bond, 



(2) 



0, otherwise, 
and we take {C<r/>} as independent, identically dis- 



tributed random variables with the distribution 

Q(0 = q6((-l) + (l-q)S(0. (3) 

The probability q that a given bond is catalytic equals 
the mean density of the catalytic bonds. The two limiting 
cases, q = and q = l, correspond to an inert substrate 
and to a homogeneous catalytic one, respectively. 

The condition of instantaneous reaction A + B — > 
is formally equivalent to allowing a NN A — B repulsive 
interaction of strength A — > oo for A — B pairs connected 
by catalytic bonds. Hence, in thermal equilibrium and for 
a given configuration ( = {C<ij>}i the partition function 
of such a two-species adsorbate is 



Z N (C) = lim Vexp[-/?H A (C)] 

A— >oo £ — 4 



(4) 



where = fc^T is the thermal energy, while the Hamil- 
tonian H\(C) = H\(() + Ho naturally separates into a 
disorder-dependent part, 



H \{Q = ^^2(<ij> {nifrij + rijirii) . 



(5) 



where the summation extends over all pairs (ij), and a 
disorder- independent contribution 

Hq = — [JAfiirij + Jsmimj + Jab (nirrij + njirii)] 
(ij) 

N 

- (/M^i + VBrrii) ■ (6) 
»=l 

In what follows we shall focus on situations in which 
the disorder in the placement of the catalytic bonds is 
annealed. In this case the thermodynamics of the system 
is given by the disorder-averaged pressure (in units of the 
lattice cell area), 



P = P(T, va,Vb) 



where 



denotes the average over all possible realiza- 



tions {C<ij>}- Once P is known all other thermodynamic 
quantities of interest can be obtained by differentiating 
P with respect to [ia, fJ-B, or T. 

Averaging Zn(() in Eq. (J2J is straightforward, 
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and yields the "effective" Hamiltonian 
H e = - ^ { [ j ab + P^ 1 ln(l - q)] (niirij + n^rrn) 



(ij) 



N 



+ JAniHj + JBiriimj} - '^{n A n l + ii B mi) ■ (9) 



1=1 



Introducing the "spin" variables o^ S {0,±1}, 



site i occupied by A (B), 




site i empty, 



(10) 



such that nj = {oi + of )/2 and rm = (— g, l + of )/2, T~t e 
can be cast into the form of the classical Hamiltonian of 
the general spin S = 1 model |l2^ . 



Tie — ~ J ^ (TiCTj — K S ' O'i&j 

<ij> <ij> 

N 

-C £ + <rrf) - H ^2 a i + A ^2 a t (11) 



iV 



1=1 



with coupling constants 



J = 
if = 



Ja + Jb~ 2 J ab _ ln(l - g) 
4 2/3 ' 

Ja + Jb + 2J AB , hr(l - <?) 



2/3 
Ma - Ms 



Ma + Ms 



(12a) 
(12b) 
(12c) 



Thus, in the case of annealed disorder, the A + B — > 
reaction model under study can be mapped exactly 
onto the GS1 model, which has been extensively analyzed 
\l2\ . The accumulated knowledge on its critical behavior, 
phase diagrams, as well as low- and high-temperature 
expansions can be straightforwardly used to elucidate 
the equilibrium properties of the present CAR model for 
general values of ma ; Ms, J a, Jb, Jab, and q, as well as 
for different types of embedding lattices 0] ■ 

In the remaining part of this paper we focus on the 
symmetric case MA = Ms and J A = Jb, implying 
C = H = so that the model reduces to the original 
Blume-Emery-Grifhths (BEG) model Q. Additionally, 
we set Jab — and consider a honeycomb lattice and a 
particular relation between K and J, for which the 2D 
BEG model, and hence the monomer- monomer reaction 
model under study, can be solved exactly |l5l fl6| . 

Following Ref. 15], in the subspace 



- pK = cosh(/3J), 



(13) 



the partition function of the 2D BEG model on the hon- 
eycomb lattice may be expressed in terms of the partition 
function of a zero-field Ising model on the honeycomb lat- 
tice, which is known in closed form In this subspace 



the 2D BEG model exhibits an Ising type phase transi- 
tion with a line of critical points obeying 



tanh(/3J) = 



2\/3 



For our CAR model Eq. [fl"3")l implies 

PJ A = (3Jb = In (1 + q) 



(14) 



(15) 



In the subspace defined by Eq. l(T5|) the disorder-averaged 
pressure (Eq. 10) is given exactly by 



(l + 2e^)--cosh((3J)-f(J), (16) 



where /i = [1a — Ms, PJ = tanh 1 (2g/(2 + e - ^)), and 
f(J) denotes the known free energy of a spin- 1/2 Ising 
model with NN interaction J on a honeycomb lattice 
(see Ref. H3)- Since in this particular case the exact 
expression for the average "magnetization" Mq = (at) is 
also known [l^J , the average densities of A and B species 
are straightforwardly calculated as (rii) = (Mq + (of ))/2 
and (mi) = (-M + (<7f})/2, where (of) = zd z (f3P) and 
z = e^ M . Furthermore, the line of critical points as func- 
tion of q (within the subspace defined by Eq. (|15|) ~) where 
a continuous transition takes place is given by 



/3mc = -M2 (qVS-1)]. 



(17) 



We now emphasize several features of these results, 
(i) For /i below its critical value (Eq. Ijl7|l ) we find 
(rii) = (mi) > (see also Fig.[2J). Upon exceeding fj, c (by 
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FIG. 2: Average density of A particles, (rii), and of B par- 
ticles, (rm), as function of their common fugacity z above 
(q = 0.60) and below (q = 0.57) the threshold value q = 1/V3 
for the concentration of catalytic segments. For clarity, the 
curves corresponding to q = 0.60 have been symmetrically 
shifted up ((rii)) and down ((rm)), respectively. The up- 
per inset (logarithmic scale) illustrates the scaling behavior 
(rii) — (m-i) ~ (z — Zc) 1 ' 8 for z — > z c . The lower inset shows 
the critical line z c (q) (Eq. 117t h the arrow indicates a path 
of increasing z at fixed q which crosses the transition line. 
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increasing the vapor pressure in the reservoirs) one of the 
densities (with equal probability) decreases sharply but 
continuously to zero while the other one rapidly attains 
unity. This reveals a spontaneous symmetry breaking 
and implies that the substrate becomes poisoned, i.e., 
most of it is covered by either one of the species. If the 
chemical potentials fiA and [Ib differ slightly, the transi- 
tion to the poisoned state is smeared out but remains de- 
tectable, (ii) The transition can occur only if the mean 
density q of catalytic bonds q is sufficiently high, such 
that q > q — 1/V3 w 0.577. For q < q, (m) = (m,) for 
all and both tend to 1/2 as (3[i — > oo. (iii) /i c > for 

q G "73 ' ' which means that in this range of q values 
the transition occurs in situations in which adsorption 
on the substrate is favored. For q e 1 , /i c < and 
hence the transition takes place for the case that desorp- 
tion into the reservoir is favored, (iv) There occur large 
scale critical fluctuations of the densities of adsorbed A 
and B particles upon approaching z c = exp(/3/x c ) from 
above or below (by varying the vapor pressure in the 
reservoirs), and the compressibility of the adsorbed phase 
diverges as \z — z c \~ 7 / 4 for z — > z c and (j,a = P>b, and as 
\fJ>A - MsT 14/15 for [i a — > and z = z c [Tg. (v) An 
analysis of the model on a Bethe lattice (coordination 
number 7 = 3) shows that the case e /3K cosh(/3 J) = 1 is 
not singular, i.e., the transition line discussed above per- 
sists for e /3K cosh(/3 J) < 1 and for 1 < e f3K cosh(/3 J) < 2 
[l3j . We finally remark that one should not expect big 
differences in the critical behavior in the case of quenched 
disorder because the transition occurs only for sufficiently 
large values of q . 

In conclusion, this study presents an exactly solvable 
model of a monomer-monomer A + B — ► reaction on a 
2D random catalytic substrate. This exact solution has 
been obtained via a mapping of the partition function 
of the two-species adsorbate onto the partition function 
of a general spin 5 = 1 model and by noticing that for 
certain relations between the corresponding coupling con- 
stants the latter reduces to an exactly solvable 2D BEG 
model 15] . In this case wc have determined the annealed 
disorder-averaged equilibrium pressure of the two-species 
adsorbate and have shown that the system under study 
exhibits a second-order (- robust in parameter space -) 
2D Ising-like phase transition if the mean density of the 
catalyst is sufficiently large. 
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